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Abstract 

Based on a recent chiral approach to nuclear matter we calculate the in-medium inter- 
action of nucleons at the Fermi surface |pi,2| = kf. The isotropic part of this quasi-particle 
interaction is characterized by four density-dependent (dimensionful) Fermi-liquid param- 
eters: fo(kf), fb(kf), go{kf) and g' (kf). In the approximation to l7r-exchange and iterated 
l7r-exchange (which as such leads already to a good nuclear matter equation of state) we 
find a spin-isospin interaction strength of g' (2m 7r ) = 1.14 fm 2 , compatible with existing 
empirical values. The consistency relations to the nuclear matter compressibility K and 
the spin/isospin asymmetry energies serve as a check on our perturbative calculation. In 
the next step we include systematically the contributions from 27r-exchange with virtual 
A(1232)-isobar excitation which have been found important for good single-particle prop- 
erties and spin-stability of nuclear matter. Without any additional short-distance terms 
(contributing proportional to A; 2 -) the spin-dependent Fermi-liquid parameters go(kfo) an d 
g'o(kfo) come out far too large. Estimates of these short-distance parameters from real- 
istic NN-potentials go in the right direction, but sizeable enhancement factors are still 
needed to reproduce the empirical values of go(kfo) and g' (kfo). This points towards the 
importance of higher order iterations subsumed in the induced interaction. We consider 
also the tensor part of the quasi-nucleon interaction at the Fermi surface. In comparison 
to the leading l7r-exchange tensor interaction we find from the 27r-exchange corrections 
almost a doubling of the isoscalar tensor strength ho(kf), whereas the isovector tensor 
strength h' (kf) is much less affected. These features are not changed by the inclusion of 
the chiral 7r./VA-dynamics. The I = 1 tensor Fermi-liquid parameters h\(kf) and h'^kf) 
follow a similar pattern. 

PACS: 12.38.Bx, 21.30.Fe, 21.65.+f, 24.10. Cn 
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1 Introduction and preparation 

In recent years a novel approach to the nuclear matter problem based on effective field theory (in 
particular chiral perturbation theory) has emerged. Its key element is a separation of long- and 
short-distance dynamics and an ordering scheme in powers of small momenta. At nuclear matter 
saturation density po — 2/^ /37r 2 ~ 0.16 fm -3 the Fermi momentum kfo and the pion mass m n 
are comparable scales (fc/o — 2mA, and therefore pions must be included as explicit degrees 
of freedom in the description of the nuclear many-body dynamics. The contributions to the 
energy per particle E(kf) of isospin-symmetric (spin-saturated) nuclear matter as they originate 
from chiral pion-nucleon dynamics have been computed up to three-loop order in refs. |2] . 
Both calculations are able to reproduce correctly the empirical saturation point of nuclear 
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matter by adjusting one single parameter (either a contact-coupling g + gi ~ 3.23 1 or a cut- 
off A ~ 0.65 GeV [2]) related to unresolved short- distance dynamics. 1 The novel mechanism 
for saturation in these approaches is a repulsive contribution to the energy per particle E(kf) 
generated by Pauli blocking in second order (iterated) one-pion exchange. As outlined in section 
2.5 of ref. [2] this mechanism becomes particularly transparent by taking the chiral limit m n = 0. 
In that case the interaction contributions to E(kf) are completely summarized by an attractive 
fcj-term and a repulsive kj-term where the parameterfree prediction for the coefficient of the 
latter is very close to the one extracted from a realistic nuclear matter equation of state. 

In a recent work |3j we have extended this chiral approach to nuclear matter by including 
systematically the effects from two-pion exchange with single and double virtual A(1232)-isobar 
excitation. The physical motivation for such an extension is threefold. First, the spin-isospin- 
3/2 A(1232)-resonance is the most prominent feature of low-energy 7rA^-scattering. Secondly, 
it is well known that the two-pion exchange between nucleons with excitation of virtual A- 
isobars generates the needed isoscalar central NN-attraction j3] which in phenomenological one- 
boson exchange models is often simulated by a fictitious scalar "a" -meson exchange. Thirdly, 
the delta-nucleon mass splitting A = 293 MeV is of the same size as the Fermi momentum 
kfo ~ 2m n at nuclear matter saturation density and therefore pions and A-isobars should both 
be treated as explicit degrees of freedom. A large variety of nuclear matter properties has been 
investigated in this extended framework in ref.jS] and it has been found that the inclusion of 
the 7rAf A-dynamics is able to remove most of the shortcomings of previous chiral calculations 
j3 El of nuclear matter. In particular, the momentum-dependence of the (real) single-particle 
potential U(p, kfo) near the Fermi surface p = kf improves significantly. The effective nucleon 
mass M* = 0.88M and thus the density of states at the Fermi surface are now better described. 
As a consequence the critical temperature of the liquid-gas phase transition is lowered to the 
realistic value T c = 15 MeV. Moreover, the isospin properties improve also substantially by the 
inclusion of the chiral 7rA^ A-dynamics, as exemplified by the density dependence of the neutron 
matter equation of state E n (k n ) and the asymmetry energy A{kf) (see Figs. 10,11 in ref. 0). 

Given the fact that both groundstate and single-particle properties can be well described by 
perturbative chiral 7rA^-dynamics up to three-loop order, it is natural to consider in a further 
step the in- medium interaction of quasi-nucleons (at the Fermi surface). Following this program 
we will therefore calculate in this work the density dependent Fermi-liquid parameters as they 
emerge from perturbative chiral pion-nucleon dynamics. 

Fermi-liquid theory was invented by Landau [6 j to describe strongly interacting (normal) 
many-fermion systems at low temperatures. At low excitation energies the elementary excita- 
tions of a many-fermion system are long-lived quasi-particles which in a certain sense interact 
weakly. The quasi-particles can thought of as free particles dressed by the interactions with 
the dense many-body medium. Landau's Fermi-liquid theory has been applied successfully to 
liquid 3 He |7] , nuclear matter and nuclei jH] . In the latter case a set of Fermi-liquid parameters, 
describing the particle-hole interaction, is assigned to nuclei heavy enough to develop a central 
region of saturated nuclear matter. Via collective excitations of nuclei and magnetic multipole 
transitions, in particular the giant Gamow- Teller resonance, some of the nuclear Fermi-liquid 
parameters [5] can be determined from experimental data. Others are related to bulk properties 
of nuclear matter, such as its compression modulus K = kj E"(kf ) or its (isospin) asymmetry 
energy A(kfo). Various calculations of the nuclear Fermi-liquid parameters using Brueckner- 
Bethe theory and phenomenological nucleon-nucleon potentials have been performed in the 

1 The cut-off scafe A serves in ref.[2] the purpose to tune the strength of an emerging attractive zero-range 
NN-contact interaction. 
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past In ref. the role of p- meson exchange for the spin-isospin parameter G' has 

been emphasized. More recently, Schwenk et al. [13.] have employed the renormalization group 
motivated universal low-momentum nucleon-nucleon potential Vi ow _k [H|. In their work the 
induced interaction generated by the particle-hole parquet diagrams plays a very important 
role. In contrast to this rather sophisticated approach, ref . ^3] takes into account only the bare 
potential V| ow _k or the in-medium G-matrix. The resulting Fermi-liquid parameters are then 
not in good agreement with empirical values (see Table 2 in ref.[15j). Moreover, saturation of 
nuclear matter could not be obtained with the potential Vi w-k in the Brueckner-Hartree-Fock 
approximation. According to recent work of Bogner et al. PB] the inclusion of the leading 
three-nucleon force from chiral effective field theory is essential in order to achieve reasonable 
saturation properties in this approach. 

The purpose of the present paper is to explore the role of the long-range two-pion exchange 
for the nuclear Fermi-liquid parameters. We will restrict ourselves to the isotropic part of the 
quasi- nucleon interaction for which empirical information is best available [H1E1E|- To be 
specific, consider the angle-averaged quasi-particle interaction V e s of two nucleons on the Fermi 
sphere, \pi\ = \p 2 \ = kf. It is irreducible in the direct particle-hole channel and it has the 
following form: 

Fo(kf) = J dVtidVt 2 (pi,P2\Ves\pi,p 2 ) 

= fo(kf) +9o(k f )a 1 ■ a 2 + fo(k f ) r x ■ f 2 + g' (k f )a 1 ■ a 2 n ■ f 2 , (1) 

where \pi,p 2 ) stands for a properly anti-symmetrized two-nucleon state and a\ >2 and T\, 2 de- 
note the usual spin and isospin operators of the two nucleons. The density-dependent functions 
fo(kf), g (kf), fb(kf) and g' (kf) of dimension fm 2 are the isotropic (/ = 0) Landau parameters. 
The spin- and isospin independent part foikf) is alternatively given by the second functional 
derivative of the energy density with respect to the occupation density in momentum space 
[HI] . Because of this property the Landau parameter fo{kf) can be directly expressed in terms 
of the angle-integrated two-body and three-body kernels, K, 2 (pi,p 2 ) and JC 3 (pi,p 2 ,p 3 ), intro- 
duced in ref.[5j to facilitate the finite temperature calculation. For two-body and three-body 
contributions to fo{kf) the following relations hold: 

f (k f f- b °^ = ^lC 2 (k f ,k f ), (2) 

/o(M (3 - body) = dpp[lC 3 (p, k f , k f ) + JCs{k f ,p, k f ) + JC 3 (kf, k f ,p)] , (3) 

where the prefactors originate from the (convenient) normalization of the kernels /C2,3 chosen 
in ref. 0. For an interpretation of eqs.(2,3) note that in a diagrammatic language the second 
functional derivative is constructed by opening two nucleon lines of a closed in-medium diagram 
which represents the groundstate energy density. The relations eqs.(2,3) serve as a good check 
on our analytical one-loop calculation of the isotropic Landau parameters fo{kf), g (kf), f' {kf) 
and g' (kf) to be presented in the next section. 

2 Diagrammatic calculation of isotropic Fermi-liquid pa- 
rameters 

In this section we present analytical formulas for the density-dependent Fermi-liquid param- 
eters fo{kf), g (kf), fo(kf) and g' (kf) as derived from tree-level and one- loop pion-exchange 
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Figure 1: In- medium scattering diagrams related to one-pion exchange, irreducible two-pion 
exchange, a two-body contact interaction, and a three-body contact interaction. 



diagrams. We give for each diagram only the final result omitting all technical details related 
to algebraic manipulations and solving elementary integrals. Diagrams can contribute via a 
direct term, where the ordering of p\ and P2 is the same in the initial and final state, and via a 
crossed term, where this ordering is reversed in the final state. In the latter case the negative 
product of the spin- and isospin exchange operators —(1 + S)(l +T)/4 has to be multiplied 
from the left. Here, we have introduced the (convenient) short hand notations, S = 0\ ■ a% and 
T = T\ • T2, which we will use frequently in the following. The spin-spin and isospin-isospin 
operators satisfy the relations S 2 = 3 — 2S and T 2 = 3 — 2T. We start now to enumerate the 
diagrammatic contributions to the angle-averaged quasi-particle interaction Foikf) defined in 
eq.(l). 

The crossed term from the (tree-level) one-pion exchange diagram in Fig. 1 reads: 



48/ 2 



(3-5)(3-T) 



1 - 



ln(l + Au 2 

4u 2 



k 2 
M 2 



(4) 



with the useful abbreviation u = kf/m n , where = 135 MeV stands for the (neutral) pion 
mass. As usual f n = 92.4 MeV denotes the weak pion decay constant and we choose the 
value qa = 1.3 of the nucleon axial- vector coupling constant in order to have a pion-nucleon 
coupling constant of g n N = 9aM/ f n = 13.2. Furthermore, M = 939 MeV denotes the (average) 
nucleon mass. The relativistic correction factor 1 — k 2 /M 2 in eq.(4) is not part of the diagram, 
but it stems from the squared M/i?-prefactor in the relativistic density of states, with E 2 = 
M 2 + k 2 . Consistency with the two-body kernel m ec L-(5) of ref. [Sj requires the inclusion 

this relativistic correction factor. 



The left diagram in Fig. 2 corresponds to iterated one-pion exchange. The direct term of 
this one-loop graph gives rise to the following (spin-independent) contribution: 



g^Mmi, 



512tt/> 2 

whereas its crossed term contributes in the form: 



(3 - 2T) 16m arctan 2u - 5 ln(l + 4m 2 



(5) 



3' ?r (4j 7r J u 



+(3-S) 



ln(l + 4m 2 ) — 2u arctan 2u + 



u arctan 2x — arctan x 

ax 

1 + 2x 2 



(6) 



Note the large scale enhancement factor M which stems from an energy denominator equal to 
the difference of small nucleon kinetic energies. The expressions in eqs.(5,6) do not include the 
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Figure 2: Iterated one-pion exchange and medium modifications of one-pion exchange and 
two-pion exchange. The short double-line symbolizes the filled Fermi sea of nucleons, i.e. the 
medium insertion —8(kf — \p\) in the in-medium nucleon propagator [2]. Effectively, the medium 
insertion sums up hole-propagation and the absence of particle-propagation below the Fermi 
surface \p\ < kf. Mirror graphs are not shown. 



contribution of a linear divergence J °° dl 1 of the momentum-space loop integral. If regularized 
by a cut-off scale A as done in ref.[2] one gets in addition the (constant) contribution: 



T {k f ) 



g\MA 
128^ 



(13T- 15-35 + 5ST) 



(7) 



which is equivalent to that of a (momentum-independent) NN-contact coupling. It turns out 
that with respect to the constraints from the empirical nuclear matter binding energy E(kfo) 
and asymmetry energy A(kf ), one adjusted cut-off A is sufficient to represent the most general 
(i.e. two-parameter) contact-interaction at this order in the small momentum expansion. At 
higher orders there are of course more parameters to describe the short- distance dynamics (see 
eqs. (22,40) below). 2 

The remaining four graphs in Fig. 2 involve medium modifications as symbolized by the short 
double-line. The second diagram in Fig. 2 together with its mirror partner can be interpreted 
as a one-pion exchange modified by the coupling of the pion to nucleon-hole states. From its 
crossed term we derive the following contribution: 



Fo{k f ) 
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Pauli blocking acts in the planar and crossed 27r-exchange box diagrams in Fig. 2 (and their 
mirror partners). The direct term from the planar box diagram contributes in the form: 
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whereas the direct term from the crossed box diagram with Pauli blocking reads: 
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(9) 
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Note that in both cases spin-dependent terms are absent. There is also the crossed term from 
the planar box with Pauli blocking. We split its contribution to JF (A;j) into a factorizable part: 
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2 At N 3 LO, corresponding to two-loop order in NN-scattering, there are in total 24 short-distance coefficients 
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and a non-factorizable part: 
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(12) 



with the auxiliary polynomial R = Au 2 + (4x 2 — l)(4y 2 — 1). These two pieces are distinguished 
by whether the (remaining) nucleon propagator in the denominator can be canceled or not by 
terms from the product of 7riV-interaction vertices in the numerator. Finally, the last diagram 
in Fig. 2 (together with three mirror partners) represents a density-dependent vertex correction 
to the one-pion exchange. The non-vanishing contribution comes from the crossed term and 
we split again into a factorizable part: 



g\Mm n 
3tt2(8^)%5' 



lu 4 + Au 2 
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(13) 



and a non-factorizable part: 
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As can be seen from the power of m n in their prefactors all contributions in eqs. (5,6,8-14) are 
of the same order in the small momentum expansion. The (bare) one-pion exchange eq.(4) and 
the contact term eq.(7) are (formally) of lowest order in the small momentum expansion. 

The p-wave (/ = 1) Landau parameter f\(kf), following fo(kf) in the Legendre-polynomial 
expansion of the spin-isospin averaged quasi-nucleon interaction, is directly related to the slope 
of the (real) single-particle potential U(p,kf) at the Fermi surface P ITU]: 



fi(kf) 



3tt 2 dU(p, k f ) 



2k 2 



dp 



(15) 



P =k f 



Its value at nuclear matter saturation density p 
mass M* via the relation: 



2/^ /37r 2 determines the effective nucleon 



M* = M 



1 - 
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2M 2 
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n -1 



/o 



3tt 2 



fi(k f o) 



(16) 



The second term —k 2 /2M 2 in the square brackets stems from the relativistic correction 
—p 4 /8M 3 to the kinetic energy. Although it is small, we keep this correction term for rea- 
sons of consistency with our earlier works E] • 

Fig. 3 shows the p-wave Landau parameter fi{kf) as a function of the Fermi momentum kf. 
The dashed curve corresponds to the approximation to In- and iterated l7r-exchange pUIUIj- 
The full curve includes in addition the higher order contributions from irreducible 27r-exchange 
with no, single and double virtual A-isobar excitation [3]. One should note that all short- 
distance parameters related to momentum-independent NN-contact interactions drop out in the 
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Figure 3: The p-wave Landau parameter fi(kf) as a function of the Fermi momentum kf. The 
dashed line corresponds to the approximation to In- and iterated l7r-exchange. The full line 
includes in addition the contributions from irreducible 27r-exchange with no, single and double 
virtual A-isobar excitation. 



quantity fi(kf). The remaining possible short-distance contribution fi(kf)^ = —hn 2 B^k 2 j/M A 
from the momentum-dependent NN-contact interactions (see eq.(9) in ref.jS.) plays numerically 
no role, since the adjustment to the saturation curve E{kf) and the potential depth {7(0, fc/o) 
gave as an optimal value = 3J. The curves in Fig. 3 display therefore the pure long- 
range effects from chiral one- and two-pion exchange. The upward bending of the dashed curve 
above kf ~ 100 MeV drives the p-wave Landau parameter fx(kf) to a large positive value of 
/i(2m 7r ) — 1.4 fm 2 at saturation density which translates into an unrealistically high effective 
nucleon mass of M* ~ 2.9M at the Fermi surface. As one can see from the full line in Fig. 3 this 
wrong trend gets with increasing Fermi momentum kf suppressed and finally reversed by the 
inclusion of the chiral 7riVA-dynamics. The negative value fi(kfo) — — 0.4fm 2 at saturation 
density kfo = 261.6 MeV [3] is compatible with the empirical values of fx{kfo), collected in 
Table 1. As a consequence the effective nucleon mass at the Fermi surface takes on now a 
more realistic value: M* = 0.88M. For Fermi momenta as large as kfo ~ Im^ the higher 
order corrections from the 27r-exchange with virtual A-excitation turn out to be essential for 
good single-particle properties. This feature which is exhibited here very clearly in Fig. 3 is in 
agreement with the findings of ref.jS]. 





k [fm 2 ] 


h [fm 2 ] 


ft [fm 2 ] 


9o [fm 2 ] g' [fm 2 ] 


ref.[Hj 


-0.3.. .0 


-0.75... 


0.9... 1.5 


small 1.6 . . . 1.7 


ref.[17j 


0.1 


-0.6 


0.7 


1.12 1.41 


ref.[13j 


-0.29 


-0.92 


0.77 


0.16 ±0.3 1.1 ±0.2 



Tab.l: Empirical values of the (dimensionful) Landau parameters at nuclear matter satu- 
ration density kfo = 263 MeV in units of fm 2 . The dimensionless Landau parameters Fq, Fx, 
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Fq, G , G' have been divided by the density of states at the Fermi surface N = 27i~ 2 kf M* 
taking for the effective nucleon mass M* = 0.8M [T2j and M* = 0.72M [THj, respectively. In 
the first row we used Nq = lfm -2 (corresponding to M* = 0.78M). 

Next, we discuss the results for the isotropic Landau parameters fo(kf), go{kf), fb{kf) and 
g' (kf) in the approximation to In- and iterated l7r-exchange. By adjusting the single cut- 
off scale A to the value A = 611 MeV this approximation leads already to a good nuclear 
matter equation of state with the saturation point at p = 0.173 fm -3 , E(kfo) = —15.3 MeV, 
a nuclear matter compressibility of K — k 2 E"(kfo) = 252 MeV j2[ [T9"] . and an (isospin) 
asymmetry energy of A(2m 7r ) = 38.9 MeV. Fig. 4 shows the corresponding Landau parameters 
fo(kf), go(kf), fb(kf) and g' (kf) (i.e. the summed contributions written in eqs.(4-14)) as a 
function of the Fermi momentum kf. The strong density dependence of fo{kf) including a 
sign change in the vicinity of saturation density kfo is a generic feature which is also shared 
by many other calculations [Oj. The other three Fermi-liquid parameters go(kf), fb(kf) and 
g' (kf) vary much less with the density p = 2k^/37c 2 (or the Fermi momentum kf). We read 
off from Fig. 4 at saturation density: fo(2m 7T ) = 0.55 fm 2 , go(2m n ) = — 1.44fm 2 , /q(2?71 w ) = 
2.01 fm 2 , and g' (2m n ) = 1.14fm 2 . The negative value of go(2m w ) = — 1.44fm 2 reflects the spin- 
instability of nuclear matter in this approximation, as discussed recently in ref. [20J. Although 
they lie outside their empirical ranges (see Table 1) the predicted values /o(&/o) = 0.55 fm 2 and 
fb(kfo) — 2.01 fm 2 do obey the Landau relations [HIE]] to the nuclear matter compressibility: 



K = 3k 2 f0 



2kf 

M* ' 7T 2 



(17) 



and to the (isospin) asymmetry energy: 



A(k f0 ) = k 



6M* 3n 



The discrepancy comes of course from the much too high effective nucleon mass M* ~ 2.9M 
in this approximation. The two spin-dependent Fermi-liquid parameters go(^/o) an d <?o(^/o) 
are connected by relations analogous to eq.(18) to the spin-asymmetry energy S(kfo) and the 
spin-isospin asymmetry energy J(fc/o) [SO]- We point out that in our calculation these (four) 
consistency relations hold with a numerical precision of one permille and better. It is also in- 
teresting to look at individual contributions. For example, the spin-isospin interaction strength 
gQ{2m 1T ) decomposes as g (2m n ) = (0.13 + 0.99 + 0.02) fm 2 = 1.14 fm 2 into contributions from 
(static) l7r-exchange eq.(4), iterated l7r-exchange eqs.(6,7), and medium modified terms. Inter- 
estingly, the medium modified terms eqs.(8,ll-14) sum up to a negligibly small net contribution. 
The important role of iterated l7r-exchange for the spin-isospin Fermi-liquid parameter g' {kf) 
has been stressed earlier in ref.jH]. The new element here is that this contribution is now 
evaluated consistently in a framework which leads to realistic nuclear matter binding and sat- 
uration. The predicted value g' (2m n ) = 1.14fm 2 of the spin-isospin Fermi-liquid parameter 
is compatible with existing empirical values (see Table 1). Nevertheless, it is clear from Fig. 3 
that higher order corrections from 27r-exchange with virtual A-isobar excitation are essential in 
order to obtain good single-particle properties. Therefore we will now turn to the contributions 
of the chiral 7rA^ A-dynamics to the isotropic Fermi-liquid parameters fo{kf), go(kf), fb{kf) and 
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Figure 4: The isotropic Landau parameters fo(kj), g (kj), f' {kf) and g' (kf) as a function of the 
Fermi momentum kf. Only contributions from single and iterated pion exchange are included. 

2.1 Two-pion exchange with virtual A-isobar excitation 

The two-pion exchange nucleon-nucleon potential is symbolized by the second diagram in Fig. 1. 
We separate regularization dependent short-distance parts from the unique long-range terms 
with the help of a twice-subtracted dispersion relation. Applying this procedure to the crossed 
term of the (irreducible) 27r-exchange diagram in Fig. 1 we obtain a contribution to J-"o(kf) of 
the form: 



where ImVc-, ImVFc*, ImVr and ImWV are the spectral functions of the isoscalar and isovector 
central and tensor NN-amplitudes, respectively. Explicit expressions of these imaginary parts 
for the contributions of the triangle diagram with single A-excitation and the box diagrams 
with single and double A-excitation can be easily constructed from the analytical formulas 
given in section 3 of ref.[4 . The fi- and /c/-dependent weighting function in eq.(19) takes care 
that at low and moderate Fermi momenta this spectral integral is dominated by low invariant 
7T7r-masses 2m 7T < fi < 1 GeV. The contributions to jF (/cy) from irreducible 27r-exchange with 
only nucleon intermediate states can also be cast into the form eq.(19). The corresponding 
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non- vanishing spectral functions read [2*T] : 



a/ '/i 2 — 4m 2 



4m2(l + 4^-5^)+/i 2 (23^ 



7l 



D + 



/i 2 — 4m 2 



(20) 



ImVr(^) 



f21^ 



vr/i(4/ 7r ) 4 

The subtraction constants which come along with the twice-subtracted dispersion relation rep- 
resentation of the 27r-exchange nucleon-nucleon potential effectively parameterize a two-body 
contact interaction (see third diagram in Fig. 1). The direct and crossed term from the most 
general two-body contact interaction up to order-p 2 give rise together to the contribution: 



7Y- 

M 2 



{35 3 + S (B 3 - 65 n , 3 ) + T (65 n , 3 - 3B 3 ) -STB 3 } 
m {B 5 + T(2B n , 5 -B 5 )+SB° + STBr}. 



57T 2 kj 



[22) 



Note that the ^/-independent term in eq.(22) involves only two independent subtraction con- 
stants. This is a consequence of the Pauli exclusion principle (or the Fierz-antisymmetric 
nature of the two-body contact interaction). The four (dimensionless) parameters B 3 = —7.99, 
Bn,3 = —0.95, B 5 = and B n5 = —3.58 have been adjusted in ref. [3] to empirical nuclear 
nuclear properties (such as the maximal binding energy per particle —E(kf ) = 16MeV and 
the isospin asymmetry energy A(kfo) = 34MeV). We note as an aside that the constant con- 
tribution eq.(7) linear in the cut-off A is of course now not counted extra since the parameters 
B 3 and B n3 collect all such possible terms. The remaining two parameters B£ and B£ T in 
eq.(22) are a priori not constrained by any groundstate properties of spin-saturated nuclear 
matter. We have also checked that the usual symmetries of the NN-interaction (isospin and 
Galilean invariance and Fierz-antisymmetry) do not imply a (linear) relation between the four 
parameters B 5 , B n>5 , B§ and B?/ . 

In order to reduce a too strongly repulsive p 2 -term in the energy per particle E{kf) a three- 
body contact interaction has been introduced in ref. jSj- The corresponding contribution to the 
quasi-particle interaction (represented by the last diagram in Fig. 1) is of the form: 



Fo(k f ) = f^(3-S-T-ST) 



(23) 



with A = 293 MeV the delta- nucleon mass splitting and the numerical parameter ( = —3/4. 
Here we have taken over from ref.jS] the parameterization of three-body contact-coupling 
strength. It is important to note here that the Pauli exclusion principle (together with isospin 
and rotational invariance) allows only for one single momentum-independent three-nucleon 
contact-coupling ~ (C#l/A/ 4 ) (NN) 3 . 

Fig. 5 shows additional pion-exchange diagrams with (single) virtual A(1232)-isobar exci- 
tation involving medium modifications. The first diagram in Fig. 5 can be interpreted as a 
one-pion exchange modified by the coupling to delta-hole states. Its crossed term gives rise to 
following the contribution: 



g\ m l u 

144vr 2 A/4 



(3-S)(3-T) 



2u 2 + 



2u 2 



1 + Au 2 



- ln(l + 4u 2 



(24) 



10 



Figure 5: Medium modifications of one-pion exchange and two-pion exchange with excitation 
of virtual A(1232)-isobars. Mirror graphs are not shown. 



with u = kf/m n , and we have already inserted the empirically well satisfied relation g^NA — 
3g w N/V2 for the 7riVA-coupling constant. The delta propagator shows up in this expression 
merely via the reciprocal mass splitting A = 293 MeV. Pauli blocking acts in the second and 
third (planar and crossed) 27r-exchange box diagrams with A-excitation. The summed contri- 
bution of their direct terms is spin- and isospin independent: 



16tt 2 A/ 4 



9 

Au 3 - 12m + 15 arctan 2u ln(l + Au 

2u 



(25) 



On the other hand the summed contribution of their crossed terms can be written in the form: 



Mkf) = -^§1 [ dx { 2X2 + Y2 + ^ ( 2X2 + 7y2 ) + ( 10X2 + 17y2 ) } > ( 26 ) 
with the two auxiliary functions: 

1 l + (u + x) 2 , „ s 

X = 2x In , 27 

2u 1 + (u - x) 2 ' v ; 



5x 2 -3m 2 -3 4x 2 + 3(1 + m 2 -x 2 ) 2 , l + (u + x. 



2 



Y = + i J -ln br. (28) 

Ax 16ux 2 1 + (u - x) 2 y ' 

Finally, the last diagram in Fig. 5 (together with three mirror partners) represents a density- 
dependent vertex correction to the one-pion exchange involving virtual A-excitation. The 
resulting contribution to J-o(kf) from the crossed term reads: 



arctan 2u 



68m 3 3 + Au 2 3 + 16m 2 + 144m 4 , 2 , 2N 

+35m — — In 2 1 + 4m 2 

3 4m 3 64m 7 v ; 
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ln(l + 4M 2 U. (29) 



9 + 30m 2 - 12m 4 + 112m ( 
+ 24m 1 

Fig. 6 shows again the isotropic Fermi-liquid parameters fo(kf), g (kf), fb(kf) and g' (kf) 
as a function of the Fermi momentum kf, including now the new contributions written down 
in eqs.( 19-29). One observes that the curve for fo{kf) has approximately the same strong 
density dependence as before (compare with Fig. 4). The value /o(fc/o) = 0.20 fm 2 at saturation 
density (kfo = 261.6 MeV [S]) lies now much closer to its empirical range (see Table 1). The 
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k f [MeV] 

Figure 6: The isotropic Landau parameters fo(kf), g (kf), f' {kf) and g' (kf) as a function 
of the Fermi momentum kf. Contributions from two-pion exchange with excitation of virtual 
A(1232)-isobars are included in addition. The (spin-dependent) short-distance parameters have 
been set to the values B§ = — 14 and B§ T = — 11. 

remaining overestimation of /o(^/o) is a consequence of the somewhat too high nuclear matter 
compressibility K = 304 MeV obtained in ref.jSj- The curves for fo{kf) in Figs. 4,6 are also 
quite instructive since they reveal that the smallness of the empirical /o(^/o) simply results 
from a zero-crossing of fo{kf) in the near vicinity of nuclear matter saturation density. The 
isospin dependent Fermi-liquid parameter /o(/c/) varies now less with density (compare Fig. 6 
with Fig. 4) and it develops even a minimum at kf ~ 290 MeV. The resulting value at saturation 
density /o(fc/o) = 1.30 fm 2 lies within the empirically allowed range (see first row in Table 1). 
We note again that the consistency relations eqs.(17,18) hold with high numerical precision 
in our calculation. However, there are also the Pauli principle sum rules (see appendix A in 
ref . (Hj ) - Due to their non-linear character these sum rules are most likely violated in the present 
perturbative approach. 

Next, we discuss the spin- dependent Fermi-liquid parameters go(kf) and g' {kf). If one 
leaves out the /^-contributions proportional to the short- distance parameters B£ and B£ T 
completely (i.e. setting Bg = Bg T = 0) one gets considerably too high values at saturation 
density, namely g (kf Q ) = 2.7 fm 2 and g' Q (kf ) = 3.0 fm 2 . Relatively large contributions come 
from the twice-subtracted 27r-exchange potential eq.(19): (2.33 5+ 1.07 ST) fm 2 , and the three- 
body term eq.(23): (0.89 5* + 0.89 ST) fm 2 . The additional short-range dynamics encoded in 
the parameters B^ and B^ T is therefore important in the spin-dependent channels. This is in 
line with ref. [12] where the role of p-meson exchange has been emphasized. In order to get 
an estimate of the short-distance parameters B£ and B£ T we bring into play the complete set 
of four-nucleon contact-couplings written down in eqs.(3,4) of ref. [22] • This set represents the 
most general short-range NN-interaction quadratic in momenta and it involves seven low-energy 
constants Ci, . . . , C 7 . After computing the spin- dependent part of Fsikf) from the direct and 
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crossed term of the corresponding contact-potential we find: 
B l = ^(-3^i + 3C 3 + 3C 4 + C 6 + C 7 ) 



30tt 2 
M 4 
320vr J 



- 6C(%) + 2C( 3 S 1 ) - 3C( 1 P 1 ) + C( 3 P ) + 3C( 3 Pi) + 5C( 3 P 2 ) , (30) 



5 5 r = ^(-3C 1 + 3C 3 + C 6 



M 4 

3(k 2 
M 4 

960tt 3 



6C(%) - 6C( 3 Si) + 9(7(^0 + C( 3 P ) + 3C( 3 Px) + 5C( 3 P 2 ) . (31) 



In that form we obtain from the entries of Table IV in ref. [22] (corresponding to various realistic 

NN-potentials) the ranges B% = -7.4 6.8 and B% T = -4.0 2.3. It is gratifying that 

these estimates lead to sizeable negative values of the short-distance parameters Pg and Bg T , 
which are actually needed in order to reduce the far too high values of go(kfo) = 2.7 fm 2 and 
g' (kfo) = 3.0 fm 2 mentioned before. The curves in Fig. 6 for the spin-dependent Fermi-liquid 
parameters go(kf) and g' (kf) have been calculated with B£ = —14 and B£ T = —11. This 
choice reproduces the empirical values g (kf ) = 0.34fm 2 and g' Q (kf ) = 1.15 fm 2 of ref. [TB"] . 
Another welcome feature is that go(kf) does now not fall below — 0.57fm 2 and therefore the 
stability condition 2n~ 2 M*kf g (kf) > —1 jUlED] is satisfied for all relevant densities. One can 
also see from Table 1 that the empirical value of <7o(&/o) has a large uncertainty. In comparison 
to the lower ends of the foregoing estimate from realistic NN-potentials sizeable enhancement 
factors of about 1.9 and 2.8 need to be applied to the short-distance parameters B% and B% T . 
Obviously, our treatment of the short-range dynamics leaves here large errors. It should also 
be mentioned that most realistic NN-potentials (with the exception of Viow-k jlH]) require 
Brueckner resummation to yield meaningful results in nuclear matter. 

The present calculation has still another shortcoming. According to earlier works EJ 
EJ EH] the induced interaction (which sums planar exchange particle-hole diagrams to infinite 
order) is an important contribution to the quasi-particle interaction in nuclear matter. Some 
of the one- loop diagrams in Figs. 2,5 do generate leading pion-exchange contributions to the 
induced interaction, but all the higher order iterations are consistently dropped here. In view 
of the good nuclear matter and single-particle properties obtained in the present perturba- 
tive framework one could argue that their effects in the spin-independent channels (/o,i(&/) 
and fo(kf)) are hidden in the adjusted values of the short-distance parameters. In the spin- 
dependent channels (g (kf) and g' (kf)) the missing higher order iterations subsumed in the 
induced interaction may be the reason for the large mismatch remaining from the estimates 
of Pg and Pg T from realistic NN-potentials. It should also be mentioned that the estimate 
P 5 ~ 7.8 from these NN-potentials is not consistent with the optimal value P5 = of ref. [3]. 



3 Tensor interaction 

In this section we investigate along the same lines as in section 2 the tensor part of the interac- 
tion of quasi-nucleons on the Fermi surface. The non-central (relative) tensor interaction has 
the following form 



p( 3 ^i ■ <?^2 ■ q ~ v\ ■ oW 2 ) ^2 [ h i( k f) + K( k f) t*i ' ^2 P(cos^), (32) 



4k f 
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with q = pi — p 2 the difference between two (incoming or out-going) nucleon momentum vectors 
on the Fermi sphere \p\\ = |_p2 j = kf. The angular dependence of the tensor interaction strength 
is represented in eq.(32) by a series in Legendre-polynomials Pz(cos$) of cost? = f>\ ■ p2- Again, 
we will restrict ourselves here to the (analytical) calculation to the leading (I = 0) Fermi- 
liquid parameters h (kf) and h' Q (kf) belonging to the tensor interaction. A first simplification 
arises from the observation that only crossed terms of pion-exchange diagrams can contribute 
to the tensor interaction. The different ordering of p\ and P2 in the initial and final state is 
a necessary condition for the spin-operators cr^ and the momentum transfer q = p\ — p>2 to 
survive in the diagrammatic amplitude. The contribution of the (bare) one-pion exchange is 
well known piT2]: 

ho(k f ) + T h' (k f ) = J?L(3 - T) ln(l + Au 2 ) (l - ^ , (33) 

with u = kf/m-x, and we have included the same relativistic correction factor as in eq.(4). The 
contribution from iterated one-pion exchange (see left diagram in Fig. 2) has the form: 

h (k f )+Th' (k f ) = ^^(3-5T)|2 + lni^-2 M arctan M 

1 f u arctan 2x — arctan x 1 

— arctan 2u + 2 / dx > , (34) 

u J l + 2x 2 J ' K ' 

and the next diagram where the pion couples to nucleon-hole states leads to the expression: 



h (k f ) + T h' Q (k f ) = ^r(3 - T) £ - 



dxx 1 



+ Ax 2 y 



2ux + [u — x ) In 



it — X 



(35) 



More involved is the evaluation of the tensor component from the planar 27r-exchange box 
diagram with Pauli blocking (see third diagram in Fig. 2). We end up with the following triple 
integral representation: 

& M *fr* K ^ /V f\.. C 



h (k f )+Th' (k f ) = |^(3 - 5T) J dx J dy f dz 



(u~ 2 + A){u~ 2 + B) 



) 1 x 2 \y + z\ (1 - x 2 ) 2 sign[(x - y){x - z)} 6{W) \ 
\2 + A + 5-4 (A + B-A)VW J ' ( } 



with the auxiliary polynomials A = 1 + x 2 — 2xy, B = 1 + x 2 — 2xz and W = (x — y) 2 (x — 
z) 2 — (1 — y 2 ){l — z 2 ). Finally, we have the last pion-exchange diagram in Fig. 2 with a density- 
dependent vertex correction. We split its contribution to the tensor Fermi liquid parameters 
into a factorizable part: 

h (k f ) +Th' (k f ) = R 9 *^™* 3 (3 - T) ln(l + Au 2 ) \8u A + Au 2 - (1 + Au 2 ) ln(l + Au 2 )] , (37) 

and a nonfactorizable part: 



g\Mm n f u ln(l + 4a; 2 ) - 4x 2 (uy/1 + Au 2 - Ax 2 + x) 2 



h (k f )+Th'(kf) = * A J (3-T) / dx ) ' — ! ^r-, \^ ■ (38) 

1 /j 0V f) 9Qn 2 fy J J VI + 4n 2 - 4a; 2 (1 + 4« 2 )(« 2 - x 2 ) 1 ; 

Fig. 7 shows the summed contributions eqs. (33-38) to the tensor Landau parameters h (kf) 
and h' (kf) as a function of the Fermi momentum kf. We have dropped the small relativistic 
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Figure 7: The Landau parameters h Q (kf) and h' (kf) related to the tensor interaction as a 
function of the Fermi momentum kf. Only contributions from single and iterated pion exchange 
are included. 

l/M 2 -correction to the l7r-exchange term eq.(33), since compared to the other contributions 
it is of higher order in the small momentum expansion. The solid and dashed line in Fig. 7 
correspond to the isoscalar h (kf) and isovector h' Q (kf) tensor interaction strength, respectively. 
Both curves start quadratically in kf and above kf ~ 100 MeV a more linear behavior takes over. 
At nuclear matter saturation density kfo = 263 MeV we obtain the values ho(kfo) = 4.78 fm 2 
and h' (kf ) = —0.83 fm 2 . These numbers are to be compared with the leading tensor interaction 
from (static) l7r-exchange which gives h (kf ) l7T = 2.68 fm 2 and h' (kfQ) l7r = — 0.89fm 2 . Quite 
surprisingly, the (parameterfree) 27r-exchange corrections lead to almost a doubling of the 
isoscalar tensor strength ho(kf ), whereas the isovector tensor strength h' (kfo) is much less 
affected. It is also instructive to see the individual contributions from the iterated l7r-exchange: 
(—0.90 + 1.50 T) fm 2 , the pion coupling to nucleon-hole states: (3.75 — 1.25 T) fm 2 , the planar 
box with Pauli blocking: (0.33 — 0.55 T) fm 2 , and the vertex correction: (—1.08 + 0.36 T) fm 2 . 
In comparison to these results previous works f2J 121] have found much smaller corrections 
to the one-pion exchange tensor interaction. The difference comes of course from our explicit 
treatment of the iterated pion-exchange (proportional to the large nucleon mass M) which is 
known to generate a sizeable tensor interaction (see eq.(32) in ref. |21j). 

Next, we are interested in the effects of the chiral 7riVA-dynamics on the tensor interaction 
at the Fermi surface. The long-range part of the corresponding 27r-exchange NN-potential 
(symbolized by the second diagram in Fig. 1) leads to the following contribution: 



where ImVr and ImlVr are the spectral functions of the isoscalar and isovector tensor NN- 



h (kf)+T h' (kf) = — / dfi /iln 



1 f 



J 2m„- 
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amplitudes [I]. The short-range pieces are again encoded in two new subtraction constants: 



ho(kf) + T ti (k f ) = i {2^2(1 - T)C(e 1 ) + 1(3 + T) [3C( 3 Pi) - 2C( 3 P ) - C( 3 P 2 



(40) 



From the entries in Table IV of ref . we obtain the average values C{e\) = — 3.95fm 4 and 
3C( 3 Pi) — 2C( 3 Po) — C( 3 P2) = — 20.5fm 4 , which will be used for the numerical evaluation of 
eq.(40). Additional contributions to the tensorial Landau parameters come from the diagrams 
in Fig. 5 involving A-excitations and medium modifications. The diagram where the pion 
couples to delta-hole states leads to the contribution: 



h (k f )+Th! (k f ) 



A f -(3-T) 



36tt 2 A/ 4 



ln(l + 4u 2 



Au 2 



1 + 4m 2 



(41) 



Somewhat more involved it the evaluation of the (second and third) 27r-exchange box diagrams 
in Fig. 5 with Pauli blocking. We find the following representation: 



h {k f )+Th' (k f ) 



9A m l fo ^/ 5 ^ 3 , 

96^A/J l >\ 3 + 



dx 



21 ft 2 
— Uu 

u \ 



2\2 ,2,2 
X ) + U + X 



L 2 

+ — (u(l + u 

u 

+ J dy 



(1 + u 2 - x 2 ) 2 + Ax 2 {x -u)(l-ux + x 2 )^ 



x 2 (l + U + X 



111 



1 + u 2 + x 2 — 2uxy 
with the logarithmic auxiliary function: 

1 . 1 + (u + x) 2 

L — — In 

Ax 1 + (u — xY 



1 + u 2 + x 2 + 2uxy 
1 + (u + x) 2 



(42) 



(43) 



The contribution of last pion-exchange diagram in Fig. 5 with a density dependent vertex cor- 
rection can be given in closed analytical form: 



h (k f )+Th' (k f ) 



4 3 

9A m l 



(3 — T) < - 64 ln(l + 4m 2 ) arctan 2u 



ln(l + 4m 2 



971^(4/^)4 

9 + 33m 2 - 192m 4 + 112m 
3^3 

2 / 2 4x 3 + 16m 2 + 144m 4 , 2/ „ 2 . . 
+ -(3 + 6m 2 -8m 4 ) + — ln 2 (l + 4M 2 )^. (44) 



u 



5M° 



Fig. 8 shows again the tensor Landau parameters ho(kf) and h' (kf) as a function of the 
Fermi momentum kf. We have now included the additional contributions eqs. (39-44) from the 
chiral 7riVA-dynamics and we have kept the relativistic /c 2 /M 2 -correction to the l7r-exchange 
term eq.(33). The solid and dashed lines in Fig. 8 for the isoscalar h (kf) and isovector h' Q (kf) 
tensor interaction strength are very similar to the previous ones in Fig. 7. At nuclear matter 
saturation density kfo = 263 MeV we find now h (kf ) = 4.93 fm 2 and h' (kf ) = — 0.69fm 2 . In 
comparison to the (static) l7r-exchange approximation this amounts again to almost a doubling 
of the tensor interaction strength in the isoscalar channel, whereas the isovector channel is 
much less affected. The short- distance term eq.(40) contributes to the total result at the 
level of about —11%. Its numerical value (—0.56 + 0.08 T) fm 2 resembles the "weak" p-meson 
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Figure 8: The Landau parameters h Q (kf) and h' (kf) related to the tensor interaction as a 
function of the Fermi momentum kf. Contributions from two-pion exchange with excitation of 
virtual A(1232)-isobars are included in addition. 

contribution of ref. ^2] (see Table 1 therein). We conclude therefore that the sizeable effects of 
iterated l7r-exchange on the tensor interaction at the Fermi surface are not altered by higher 
order corrections from (irreducible) 27r-exchange with virtual A-isobar excitation and explicit 
short- distance contributions. 

It is well known that the expansion in Legendre-polynomials eq.(32) converges slowly in 
case of the tensor interaction ^2j E3]. For that reason we have calculated also the I = 1 
tensor Fermi-liquid parameters h\{kf) and h[(kf) in the present framework. The corresponding 
(numerical) results are shown as a function of the Fermi- momentum kf in Fig. 9. The dashed 
lines correspond to the approximation to single and iterated pion-exchange and the full lines 
include in addition the (higher order) contributions from 27r-exchange with virtual A(1232)- 
isobar excitation. In comparison the leading contribution from (static) l7r-exchange: 



one finds a pattern similar to the I = case. At nuclear matter saturation density kfo = 
263 MeV one has now from static l7r-exchange: hi(kfo) l7T + T h[(kfo) ln = (3.33 — 1.11 T) fm 2 . 
These numbers get enhanced to hx(kfo) +T h[(kfo) = (5.08 — 1.83 T) fm 2 by the iterated pion- 
exchange effects, whereas the further inclusion of contributions from 27r-exchange with virtual 
A(1232)-excitation leads only in minor changes: hi(kfo) + Th'^kfo) = (4.37 — 1.85 T) fm 2 . 
One should also note that to the order in the small momentum expansion we are working 
here, h\{kf) and h[(kf) receive no contribution from short-distance contact terms. We can 
therefore conclude that the features observed for the / = tensor Fermi-liquid parameters are 
approximately repeated at / = 1 and presumably they continue to hold at higher /. 



h^kfY^ + TKikf) 




(3 



T) (2 + m~ 2 ) ln(l + 4u 2 ) -4 , 



(45) 
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Figure 9: The Landau parameters h\(kf) and h'^kf) related to the tensor interaction as a 
function of the Fermi momentum kf. The dashed lines show contributions from single and 
iterated pion-exchange only. The full lines include also contributions from two-pion exchange 
with virtual A(1232)-isobar excitation. 

4 Summary and conclusions 

We have explored in this work the role of the long-range 27r-exchange for the in-medium inter- 
action of quasi-nucleons at the Fermi surface \p\^\ — kf. Our analytical calculation is rooted 
in a recent chiral approach to nuclear matter which can successfully explain binding and sat- 
uration of nuclear matter through two-pion exchange mechanisms 0] • The isotropic part of 
the quasi-nucleon interaction is characterized by four density-dependent (dimensionful) Fermi- 
liquid parameters fo(kf), fb(kf), go(kf) and g' {kf). In the approximation to l7r-exchange and 
iterated l7r-exchange we find a spin-isospin interaction strength of g' (2m n ) = 1.14 fm 2 , compat- 
ible with empirical values. The consistency relations to the nuclear matter compressibility K 
and the spin/isospin asymmetry energies are fulfilled in our perturbative calculation. However, 
the non-linear Pauli principle sum rules [HJ Ej are most likely violated. 

In the next step we have included in the quasi-particle interaction the contributions from 
27r-exchange with virtual A(1232)-isobar excitation. This extension to higher orders in the 
small momentum expansion is necessary in order to obtain good single-particle properties and 
to guarantee spin-stability of nuclear matter [2*0] . Leaving out the short- distance terms (con- 
tributing proportional to kj) the spin-dependent Landau parameters g Q (kf Q ) and g' (kf ) come 
out far too large. Estimates of these short-distance parameters from realistic NN-potentials 
go in the right direction, but sizeable enhancement factors are still needed to reproduce the 
empirical values of go{kfo) — 0.34 fm 2 and g' (kf ) ~ 1.15 fm 2 |13| . Another shortcoming of the 
present calculation is that only leading one-loop pion-exchange contributions to the induced 
interaction [TQl H2J UH1 E3] could be included consistently. This is presumably the reason for 
our failure to describe the spin-dependent Fermi-liquid parameters, which are not constrained 
by any bulk properties of nuclear matter. Actually, one conclusion of ref. [2H] has been that an 
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accurate description of the spin-response should include the induced interaction. 

We have also considered the tensor part 3 Bx-qBz-q—avVi q 2 of the quasi-nucleon interaction 
at the Fermi surface. In comparison to the leading l7r-exchange tensor interaction we have found 
from the iterated l7r-exchange corrections almost a doubling of the isoscalar tensor strength 
h (kf), whereas the isovector tensor strength h' (kf) was much less affected. This feature did 
not change qualitatively by the inclusion of the chiral 7riVA-dynamics. A similar pattern has 
been observed for the I = 1 tensor ial Fermi-liquid parameters h\{kf) and h[(kf). 
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